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1 Introduction 

Scattering amplitudes are very important in the context of black hole physics. They 
relate directly to astrophysical problems but also to other theoretical problems like 
stability of gravitational solutions and AdS/CFT duality. Rotating black holes are 
particularly difficult to study and are usually treated using semi-analytical methods 
and matched asymptotic expansions. 

Recent work have pointed to the possibility of extracting exact scattering ampli¬ 
tudes (in some particular cases) using only monodromy data of the radial part from the 
wave equation of interest [1, 2]. In the latter work, the authors outlined a procedure to 
accomplish this program for conformally coupled scalars in a generic Kerr-NUT-(A)dS 
black hole. One of the results showed that the scattering amplitudes are completely 
determined by the composite monodromy parameter a, and that the theory of isomon¬ 
odromic flows [3-5] exposes a hidden, non-linear symmetry of the parameters which can 
be used to relate the problem of finding a to the connection problem of the Painleve 
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VI transcendent. In this paper, we obtain explicit analytic expressions for a and the 
scattering coefficients using recent results for the Painleve VI r-function expansion in 
terms of c = 1 conformal blocks [6, 7]. We also show that both the scattering problem 
for the radial equation and the eigenvalue problem for angular equation can be solved 
by the isomonodromy method. Finally, we obtain the hrst corrections for a in both 
near-extremal Kerr-dS limits. 

One important aspect of the master perturbation equation of spin s helds for the 
Kerr-dS background - also called Teukolsky master equation (TME) - is that it is 
not only separable into radial and angular parts, but it can be reduced to a Fuchsian 
differential equation with 4 singular points [8]. This equation has been studied since 
the late 19th century and it is called Heun’s differential equation when written in the 
canonical form: 


y' + 



6 

z-l 



y + 


a/3z — q 


z{z — l){z — t) 


y = o, 


( 1 . 1 ) 


with its coefficients obeying the Fuchs condition '-^ + 5 + e = a+(3 + 1. From this starting 
point, series expansions for scattering amplitudes in the Kerr-dS background have been 
obtained using a different method than ours by Suzuki, Takasugi and Umetsu in a series 
of papers [9, 10]. Inspired by earlier works of Erdelyi and Schmidt (for references, see 
the review book [11]), Suzuki et al. used a hypergeometric series expansion for Heun 
functions: 


OO 

yy{z)= ^ a];2Fi(-i/-n-i I + f;7,^), 

n=—oo 


( 1 . 2 ) 


to study scattering in the Kerr-dS background. This series expansion converges with 
the correct local behaviour near z = 0 and z = 1 for special values of the coefficient z/ 
coming from an augmented convergence condition. 

The interpretation of z/ is not so clear in the literature, but in the approach pre¬ 
sented here its interpretation becomes quite natural: it is associated to the composite 
monodromy (jjj of the full solution around two singular points corresponding to the 
two horizons involved in the scattering. In particular, the series solution (1.2) has the 
correct local behaviour near z = 0 - the outer horizon r = r+ - and z = 1 ~ the 
cosmological horizon r = rc- The monodromy coefficient of (1.2) at 2 ; = oo is equal to 
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the monodromy at infinity of the family of hypergeometric functions in ( 1 . 2 ), which is 
6^00 = 2z/ + l (mod2n). With respect to the hypergeometric functions, this is equivalent 
to the composite monodromy of a loop enclosing z = 0 and z = 1 simultaneously and, 
therefore, u parametrizes the composite monodromy between 0 and 1 of the full Heun 
solution ( 1 . 2 ). 

Along with the obvious relevant applications for the scattering theory of black 
holes, the astrophysical and stability studies that issue from it, the extracting of the 
connection coefficients for the Heun equation described in this article is a century-old 
problem in mathematics, deeply tied to the Riemann-Hilbert problem. In its origi¬ 
nal form, this problem posed the question of writing an ordinary differential equation 
with prescribed monodromy data. As will be explained in Section 2 , we will be pri¬ 
marily interested in the reverse Riemann-Hilbert problem, where one is interested in 
extracting the monodromy data from the differential equation. The major theoretical 
breakthroughs for solving this problem were achieved by Schlesinger^, who discovered 
the non-linear symmetry of the monodromy data, encoded in the Schlesinger equa¬ 
tions, which we will revise in Subsection 3.1. Another milestone came about with the 
work of Miwa, Jimbo and various collaborators, building up from critical systems in 
two-dimensional statistical mechanics [3-5], where the Hamiltonian structure of the 
Schlesinger equations (see Section 3) were explored to finally prove the Painleve prop¬ 
erty of the solutions [13], fostering tremendous developments in integrable systems and 
the search for similar structures in other areas. 

Another great leap came in from the AGT conjecture [14], its subsequent proof 
[15] and the long series of applications to Painleve transcendents, specially [ 6 , 7]. In 
these combinatorial solutions for the Painleve VI r-function were given exploring the 
relationship between accessory parameters of the Heun equation and four point func¬ 
tions in conformal field theories, via the Virasoro conformal blocks. We will give a 
very sketchy description of this relationship in Section 4. We will make the claim that 
this achieves the analytical solution of the scattering problem: an implicit solution for 
the scattering coefficients will be given in terms of the Painleve VI r-function and the 
radial equation parameters. Likewise, the eigenvalue problem for the angular equation 
can also be cast in terms of monodromy data, yielding a formal solution, as explained 

^Some of the historical development of the Riemann-Hilbert problem and its relation to the theory 
of Painleve transcendents can be found in [12]. 
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in Appendix B. Both results are analytical and valid for generic ranges of black holes 
and scalar wave parameters, and we list expansions for the Painleve VI r-function and 
scattering coefficients in Appendix A. The exploration of these solutions to extract 
physical intuition about black hole physics is a very interesting future problem. 

The relationship between Fuchsian equations and the theory of special functions 
in the complex plane was clear since the beginning of the former held. The obvious 
SL(2, C) symmetry can be used to reduce some of the parameters, much in the way 
one can reduce the general Riemann differential equation to Gauss’s hypergeometric 
canonical form. What is perhaps less obvious is that one can use deep results from 
the Virasoro algebra representation theory to essentially solve the monodromy and 
connection problem for generic linear systems. We discuss the implications of this for 
scattering and black hole physics in Section 5. It is worth noting that, in a parallel 
result [16], the authors gave analytical results for the standard asymptotically hat Kerr 
black hole in terms of the Painleve V r-function, itself related to irregular conformal 
blocks. 

Generic spin perturbations, described by the Teukolsky master equation [8], reduces 
to a diherential equation of the same nature as eq. (1.1) and should be amenable to 
the same methods outlined here. In fact, this simplihcation happens for every vacuum 
type-D metric [17]. Scalar helds in higher dimensional black holes backgrounds yield a 
more complicated isomonodronic structure, possibly tied to Wn conformal blocks [18]. 
More dimensions means more singular points and thus a higher-point isomonodromic 
how. This case is not so well studied as the one with 4 singular points but it is 
known that the isomonodromic how still has the Painleve property [12] and thus the 
asymptotics should be similar as in our case. Therefore, although we shall focus on 
scalar perturbations of 4-dimensional Kerr-dS black holes, our method should be useful 
to understand the same integrable structure in more general cases for higher-spins and 
higher- dimensions. 

2 Kerr-de Sitter Wave Equation 

The metric of a 4-dimensional rotating black hole with de Sitter (dS) asymptotics 
has been obtained by Garter in the late 1960’s [19] and the generalization for higher- 
dimensions, with the addition of a NUT charge, was found in [20] . The special property 
that guided Garter to hnd this metric was the separability of the Hamilton-Jacobi 
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equation for the geodesic motion. We now know that not only Hamilton-Jacobi but 
linear perturbation equations for this metric are also separable for spin 0, 1, | and 

2 [8, 21, 22], in the formalism of the Teukolsky master equation, which is true for 
any Petrov type-D spacetime. For convenience, the Kerr-dS metric can be written in 
Chambers-Moss coordinates as: 


ds^ = — 
where 
A,= 


_ a Odcf)') ^ [adt - (r^ + a^)d(j)Y + ^ 

P X P X \Ae ArJ 

= 1 + with the dS radius = 3/A, and 

1 + cos^ 9, Ar = (r^ + a^)(l — — 2Mr, cos^ 9. 


( 2 . 1 ) 


( 2 . 2 ) 


We have chosen the dS radius, but one can write the metric in terms of the AdS radius 
just by Wick rotating the radius L ^ iL and the Kerr metric can be recovered by 
making L —)■ cx). The coordinate singularities of this metric are now given by the root 
of a 4th order polynomial = 0. For a certain range of black hole parameters, we 

can hnd 4 real roots in the dS case, which we call (r_,r_,r+,rc), and 2 real roots 

in the AdS case, called (r_, r+, C, C)- The roots r_ and r+ are the inner and outer 
horizon, respectively, like in the pure Kerr case. In the dS case, one of the roots is 

usually addressed as non-physical, r _= — (r+ -1- r_ rc) is a negative number, and 

rc is the cosmological event horizon [23]. To clarify the causal structure, we show the 
Penrose diagram of Kerr-dS for 0 = 0 in Fig. 1. This diagram continues indehnitely in 
all directions. The dashed vertical line again represents the black hole singularity [24]. 

We dehne two Killing vectors: 


= dt + Q{r+)d(j), = dt + Q{rc)d(j), 

such that they are null at each respective horizon r_|_ and rc- This entails to the con¬ 
stants Qk = f^(’"fc) being the angular velocities of each horizon. In particular, this 
induces a frame-dragging effect near the event horizon, as no observer can stay station¬ 
ary with respect to dt and is forced to co-rotate with the horizon. The angular velocity 
and temperatures of the event and cosmological horizons for an observer following ^k 
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r = oo 


r = 00 



Figure 1. Causal diagram of maximally extended Kerr-dS black hole for 0 = 0. 


orbits, with k = C, are given by: 




rl + a^'’ 


T, 


\Kirk)\ 

4:7rx^{rl + a^) ’ 


(2.3) 


in which we choose the sign of to both be positive temperatures. 

Let be a solution of the Klein-Gordon equation 

for D = 4 Kerr-dS metric in Chambers-Moss coordinates (2.1). The radial equation 
resulting from this solution is: 


dr{ArdrR{r)) + 


12^ 2 X^(a;(r2 + a2) 

L2"' A, 



R{r)=CeR{r), 


(2.4) 


where the separation constant between the angular and radial equations is Ci = Xe + 
X^(a^a;^ — 2mau}). The parameter ^ is the coupling constant between the scalar held 
and the Ricci scalar. Typical values of the parameter ^ are minimal coupling ^ = 0 and 
conformal coupling .^ = 1/6. In the latter, (2.4) is equivalent to the Teukolsky master 
equation for a spin zero perturbation [8, 9]. The angular equation has essentially the 
same form as the radial one and we are also able to obtain formally the Xe eigenvalues 
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with the method outlined below, as described in the Appendix B. 

If we restrict to the conformally coupled case, the equation (2.4) can be reduced 
to a Heun equation [2, 8, 17]: 


y"{z) + 


1-00 , 1-01 , 1 - 0 ^ 
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z-l 
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^0 


z-to 


y\z) + 


m2 


z{z-l) 


where y{z) is given in terms of R{r) by: 


^o(^o — l)A'o \ 
z{z - l){z-to)) 


y{z) = 0, 
(2.5) 


y{z) = z^^l\z - \)^^l\z - Rt^l\z - z^)-^R{r)- (2.6) 


and 


Z ZoyQ , ZoyQ 

r — r _ 


rc — r _ 

rc - r- 


r_|_ — r_ 

^0 ^oo" 

r_|_ — r _ 


(2.7) 


We note that the scalar Teukolsky master equation is equivalent to the Klein-Gordon 
equation for a conformally coupled massless scalar held. The monodromy coefficients 
are 


0fc = 


uj{rl + a^) — am 

Kin) 


i / oj — Vtkm\ 

^ V n ; ’ 


/c = 0,1, to, C)0, 


( 2 . 8 ) 


where we use plus or minus sign to have positive temperatures T^, and 


Kn = 


to 


1 + 


r+ 


r_ 


-2fx 


A;(r+) 

2 Ci;(r_|_r_+ a?) 


—jri + \e + — 2aum) 


am 


A;(r+) 


(2.9) 


The values of 0^ obey Fuchs relation 0o + 0i + Oto + Qi + <12 = and q 2 — Qi = 0oo- 
Also, in terms of (2.5), we have that qiq 2 = 1 + 0oo- The set of 7 parameters 
(00, 01, 0*0, ^ 1 ,A'o) in (2.5) are related by the Fuchs relation, and we see that 
the resulting 6 parameters dehne the Heun equation and its fundamental solutions. 
The Heun equation has a rich history in mathematical physics. For such details about 
it, we refer to [11, 25]. 
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The Riemann symbol for (2.5) is 


0 1 to oo 

0 0 0 gi z 

^0 (^to Q2 


yi~) = P 


( 2 . 10 ) 


and the parameters (2.8) and (2.9) are complex, so it is not trivial that, given a solution 
y{z), the complex conjugate {y{z))* will also be a solution. However, it can be checked 
that {y{z))* satishes (2.5) with parameters (— 6 * 0 , —Oi, —OtQ ; to,- 6 ^ 0 )) when the physical 
parameters are real - note that to is real in the de Sitter case. More importantly, the 
radial part of the held R{r) behaves in a well-determined manner under time-reversal 

= R-ujA-m{r)- ( 2 . 11 ) 

Thus, y{z) and its complex conjugate are also related by time-reversal. This can be 
checked by inspecting the transformation that brings the radial equation (2.4) to the 
canonical form (2.5). If one invert the signs of the in (2.6), one arrives at the 
ODE satished by {y{z))*. Hence, time-reversion symmetry tells us that a set of linearly 
independent solutions of (2.5) is: 

y{z) and z^°{z-lf^{z-toY^o[y[z)y. ( 2 . 12 ) 

Note that the Riemann symbol for both solutions are the same. This fact will be im¬ 
portant below when trying to compute the scattering coefficients in terms of connection 
data. 

3 Isomonodromic Approach to the Scattering Problem 

Consider now two linearly independent solutions y^^\z) and y^‘^\z) of (2.5). The con¬ 
nection problem for the ODE consists in writing a solution with known behavior near 
one singular point, yl^\z) at z = Zi, as the particular linear combination of solutions 
with known behavior at another critical point z = Zj, 


!/th) = (Eij)nyf\z) + {Etj)i2yf{z). 


( 3 . 1 ) 


Along with the connection coefficients (-Eij)ii and (-Eij)i 2 for there are the 

supplementary connection coefficients {Eij) 2 i and {Eij )22 for the linearly independent 
solution yi‘^\z) with known behavior ai z = Zi. So the Eij should really be thought of 
as matrices. 

Solving the connection problem is intimately related to the monodromy problem. 
For Fuchsian equations, the independent solutions with known behavior at z = Zi are 
of the form: 

= (J - + 0{z - z,)) (3.2) 

■fl 2| 

where a\ ’ ^ are the solutions of the indicial equation at Zi. Because in general these 
are non-integers, we have that — Zi) Zi) = ^{z). The matrix 


g2^iai Q \ 

0 ) 


(3.3) 


implements the effect of this monodromy around the solutions with known behavior at 
Zi- If, however, we decide to write the matrix Di not in terms of {z) but in terms 

f 1 2^ 

of 2 /j ’ ^{z), we will have to conjugate it using the connection matrix: 

Mij = E,jDj{Eij)-\ (3.4) 


where Mij is now the matrix associated to the monodromy around Zi written in the 
natural basis for Zj - he., those obtained by the Frobenius method from the solutions 
of the indicial equation, where the monodromy around Zj is diagonal. 

It is clear now that if we use another basis, not necessarily one associated to a 
critical point, the matrix Mj associated to the monodromy around Zi will be related 
to Mij by conjugation. Solving the monodromy problem amounts to hnd the set of 
Mi associated to monodromies around all singular points of a known ODE. This is 
the reverse of the Riemann-Hilbert problem (in its initial guise), which is to write an 
ODE with a known set of monodromy matrices Mj. For Fuchsian equations with up 
to three regular singular points, the Mj - and hence the Eij up to normalization - are 
completely determined from the solutions of the indicial parameters Qi , information 
readily available from the ODE. For four or more regular singular points, the Mi depend 
in a very non-trivial way on the accessory parameters, which in the Heun equation case 
(2.5) is essentially related to to and Kq. When there are irregular singular points, Stokes 
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parameters also come into play [1, 3]. This happens for the scalar wave equation of the 
Kerr black hole and we treat this problem in a separate paper [16]. 

Also, a simple parameter counting argument shows that we cannot have all pos¬ 
sible sets of Mi represented by the initial ODE - we have to introduce apparent sin¬ 
gularities [12], So, instead of studying the Riemann-Hilbert problem from the ODE 
perspective, it is more illuminating to use a matricial system: 

= A{z)Hz), (3.5) 

where 



with Ai being matrices independent of z. Each row of the matrix 


/ 1 /( 2 ) ( 2 ;)\ 


(3.7) 


consists of linearly independent solutions of the ODEs below. This matricial system is 
generically called a Fuchsian system. It can be easily verihed that any element of the 
hrst row of $( 2 :) will satisfy the ODE: 


dly - (Tr A -h log Ai 2 )d^y -1- (det A - ^^An -F A^d^ log Ais)// = 0, (3.8) 


whereas elements of the second row will satisfy a similar equation, but with indices 1 
and 2 interchanged: 

dlw - (Tr A -F dz log A 2 i)dzW -t- (det A - i9^A22 + A 22 dz log A 2 i)tc = 0. (3.9) 

It is also straightforward that the elements of a given row will be linearly independent. 
The respective Wronskians for (3.8) and (3.9) are: 

Wi{z) = Au(z) Det $(^), W 2 {z) = A 2 i(z) Det <I)(; 2 ), (3.10) 

and theorems of existence and unicity of solutions tells us that any two solutions <I>(^) (z) 
and <I>(2)(2;) of the Fuchsian system are related by right multiplication by a constant 
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matrix g, i.e, 

The Fuchsian system connects the Riemann-Hilbert problem with the theory of 
flat holomorphic connections [26]. From equation (3.5) for ‘F(^), we see that: 




( 3 . 11 ) 


which suggests the interpretation of A{z) as a flat GL(2,C) connection. Clearly, ‘F(z) 
is uniquely dehned by an initial condition <F(zo) = *^*0 and changing the condition 
amounts to a conjugation transformation $( 2 :) = ^{z)g, for g G GL(2,C). The formal 
solution for <F(z) above is 


$( 2 :) = Pexp ( / A{z)dz] $ 0 , 


’zo 


( 3 . 12 ) 


therefore, the poles of A{z) correspond to the branch points of $( 2 ;). We thus see that 
the monodromy problem can be recast as a holonomy problem of the connection A{z). 

We will now set the parameters of Ai so that the equation (3.8) mimics the Heun 
equation (2.5). Given the system (3.5) with 4 regular singular points in the canonical 
gauge, we parametrize the Ai following [4]: 


A=l i = 0 .M. 

,5(p.+ »<) -Pi 


SO that 


Aoo — “(^do -l- Ri -l- At) — 



ei = TTAi, et = TTAi, 


and 


keC, 


( 3 . 13 ) 


( 3 . 14 ) 


( 3 . 16 ) 


2;(2; — 1)(2; — f) 

with A being a complicated function of the Pi and g*. With this parametrization, we 
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can find that each hrst row element of *h(^) satisfy 


y" +p{z,t)y' + q{z,t)y = 

, , 1-00 1-01 1 - 0 i 1 

p{z,t) = -^ +^ + 


z — 1 z — t z — A’ 

ki{k 2 + 1) t{t-l)K A(A-l)/i 

q[z,t) = — --^^ + 


z{z — 1 ) z{z — l){z — t) z{z — l){z — X)^ 


(3.16a) 

(3.16b) 

(3.16c) 


where p is Aii[z) calculated at 2 ; = A. By setting A = t, 0t = 0*^ — 1 and —p + K = Kq, 
we can recover (2.5) by making t = to- An explicit parametrization is given by 

(00 + 01 — 0t + 0oo)(0O — 01 + (1 + 2to)0i + 0oo) to(to — 1) f 000* 010* 

Po = - 7 ;;-^- 7 - 1 Ao + ^-h 


40. 


0 . 


io ^0 — 1 / ’ 


(00 + 01 — 0* + 0oo)(0O — 01 + (1 + 2to)0* — 0oo) to(to — 1) 000* 010* 

Pi = --7- Ao + --h 


(At)21 — 


4000 

Po + 00 Pi + 01 


0O 


^0 ^0 — 1 / 


Pi 


Qo 


qi 


Pt = -0*, Qo = -01, 0* = 0, 

Po 


(3.17) 

where 6^0 = Ki — K 2 - Even with these choices, there is still some freedom in choosing 
the elements of A{z) given the ODE, since the value of go and qi is determined up to 
a multiplicative factor. From the ODE (3.8), it is clear that 


A{z) and e-'''"M(7)e'''"^ 


( 3 . 18 ) 


with q some arbitrary complex constant, will yield the same ODE for both rows of <h. 
In terms of the parametrization (3.13) this corresponds to an overall scaling of the gds, 
or still a rescaling of k in ^412. We will use this freedom below to £x the normalization 
of the solution. Note that this transformation maintains the Wronskian invariant. 

The boundary conditions associated with the Fuchsian system with the choice 
(3.14) are schematically given by 


iGi{l +0{z - Zi)){z - Zi)^0 0 )g-^ Z ^ Zi, 

Hz) = { 0 ^ 

I (1 + 0{z ^))z V 0 ^ 2 /^ 2 ; —)■ CX), 


(3.19) 


- 12 - 



















where 


= P,20) 

\ qiPi / 

and Qi are the connection matrices, to be discussed below. Explicitly, we have Gi as 
the most general matrix such that = Gj ( q) With this parametrization, we 
dehne the canonical - or “natural” - solutions near Zi to be ^i{z) = ^{z)g~^ and note 
the leading behavior for each entry of ^i{z), based on the parametrization for Gi, 


yi^\z) = Xie^^qi{z - Zif^{l + ...), 
wl^\z) = Xie^*{z - ZiY^il + ...), 


= Xie 
wi^\z) = XiC 




qiPi 


(3.21) 


where we are disregarding 0{z—Zi) subleading terms in both Frobenius expansions. The 
values of Aj are determined from boundary conditions once we £x that the connection 
matrices should satisfy det gi = 1. Given that one can choose boundary conditions such 
that 

W{z) = det = z^'^iz - (3.22) 

then, for det gi = 1, we compare the asymptotic expressions (3.21) and hnd 


'-'i 


(3.23) 


/1 o\ 

We will £x the real part of the /?* by requiring that yl ’ {z) in (3.21) are associated 
with equal but opposite radiation fluxes. The radial equation has real coefficients, then 
if R{r) is a solution of (2.4), so is its complex conjugate R*{r). One can then easily 
hnd that the quantity 

^= -'Mr mm - MM 

is independent of r if R{r) is a solution of (2.4). Physically, this corresponds to the 
radiation hux in the r direction, with a prefactor chosen for later convenience. In terms 
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of y{z) and y*{z), (3.24) becomes 


J = if{z) 


W{zf/S\z)j-^{W{z)-^/^y{z)) 

-W{z)-^/^y{z)j-^{W{zY/^y*{z)) 


W' 


= if(z) { W[y,y*] - ’-^\y\^ 


(3.25) 


(3.26) 


where f{z) = z{z — 1)(2; — to) and W[yi,y 2 ] = yiy '2 — y 2 y'i is the Wronskian of two 
functions. Near 2 ; = 1, we hnd that 


p/3i+/3r _ ^ 

kil’ 


(3.27) 


and thus the value of 77*-^^ = —is 


77 ( 2 ) = \p^q^ 


(3.28) 


which is independent of z and then will allow us to compare the normalization of the 
solutions yf''‘^\z) at different points. One can repeat this calculation for the solutions 
near z = t, for example, and check that the same answer as above is obtained. 

The scattering problem is formulated in terms of the normalized radial wavefunc- 
tions uf, for example, 

M7(r) = ^u^{r) + p<(r), (3.29) 

where uj{r) represents a normalized purely incoming wave at the black hole outer 
horizon z = t and wf (r) represent normalized incoming and outgoing waves at the 
cosmological horizon z = 1. The problem is complemented by its time-reversed version 

u+ (r) = (r) + yut (r), (3.30) 

where a purely outgoing wave at the black hole horizon divides into a superposition of 
incoming and outgoing waves at the cosmological horizon. We can now £x k in (3.15) 
and the scattering coefficients in such a way that yl^'^ corresponds to m7 {r) and yl‘^^ to 
uf{r). With this provision, we can see that the normalized connection matrix between 
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z = t and z = 1 is of the form: 


B« = a9r' = (3.31) 


which simplihes our calculations somewhat because we will not need the full set of mon- 
odromy matrices to compute the scattering coefficients. Now, consider the composite 
monodromy au dehned by 

2 cos TTCTii = Tr(MiMt). (3.32) 

The monodromy matrix around z = Zi can be written as Mi = g~^Digi. Plugging this 
into (3.32) and, using (3.31), we hnd that 


TV 


2 sin 7 i9i sin irOf 
COS7r(6*i — 6t) — COSTTCTit 


sin TiOi sin nOt 

sin |(crii + - 9t) sin |((Tu -9i + 9t)' 


(3.33) 


Therefore, the scattering amplitudes only depend on the monodromy data. When an 
is either real or purely imaginary, (3.33) is real and we have a well-dehned scattering 
problem with uf being the complex conjugate of . However, when an is complex, 
we notice that the amplitudes are dehned up to a phase. Therefore, we set 


T = 


(3.34) 


and this extra phase can be absorbed in the imaginary part of f3t in (3.21), for example, 
which is not hxed a priori by our radiation hux argument. Finally, we end with 




sin Ti9i sin 

sin |(crit + 9i - 9t) sin |((Tit -9i + 9t) 


(3.35) 


Now, the only non-trivial global information we cannot read directly from the ODE 
is the composite monodromy parameter ait, which is not easy to calculate from (2.5). 
This is the subject of the following section. 

Even without a deep study of the parameter an, one can learn some lessons from 
the general structure of the formula for the transmission coefficient (3.35). First and 
foremost, superradiance occurs for frequencies u and azimuthal angular momentum m 
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in the range where both 6 ^^^ and 9i are (imaginary) negative: 


Gta = - 

° 271 T+ 


i u) — 


i oj — Vtcm 


(3.36) 


27r Tc 


and quasi-normal modes are found from the vanishing of the denominator, which poses 
a quantization condition for exit, 


ait i (6^1 — 6t) = 27m, n G Z. 


(3.37) 


Similar conditions also hold for the inner-outer horizon quasi-normal modes. The im¬ 
pact of these results for the general problem of inner horizon instabilities will be left 
for future work. 

3.1 Isomonodromic Deformations 

The isomonodromic deformation of (3.8) is a non-linear symmetry with an interesting 
interpretation in terms of flat holomorphic connections and transcendental functions. 
We have previously pointed out [2] the importance of this symmetry to the analytical 
solution of the scattering problem, and refer to it for a review on the subject. Following 
[3-5, 12], if we take the partial fraction expansion of A{z) then the system 



(3.38) 


represents the components of a flat holomorphic connection in the two-dimensional 
complex space ; 2 , t, satisfying d^A — dtB — [A, 5] = 0 if A^ satishes the Schlesinger 
equations: 



(3.39) 


Given that the connection is flat, we have that, along the solutions of the Schlesinger 
equations with respect to the flow of t, the monodromy data is the same - the matrices 
Mi are mantained up to overall conjugation. This warrants the name “isomonodromic 
deformation” for the solutions of the Schlesinger equations. 

Let us write the Schlesinger equations for the corresponding differential equation 
(3.16). The singularity at z = A in (3.16) is known as an apparent singularity, because 
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its monodromy is trivial^. This implies in an algebraic constraint between K, X and 
/i, which will be explicited below. The parameters /i and A can be seen as canonically 
conjugate variables if the Schlesinger equations are written as the Hamiltonian system: 




(3.40) 


with the canonical Poisson bracket {/,S'} = 


This Hamiltonian system can be used in principle to calculate the monodromy 


data, but as it turns out (3.40) can only be solved in general in terms of the Painleve 


VI transcendent, whose general properties are still unknown [27], but are under active 
study given its relation to the AGT relation in Liouville held theory [7, 28, 29]. The 
problem immediately relevant to us is to extract the monodromy parameter au (3.32) 
from the asymptotics of the Painleve system. The system is determined by (3.40) with 
initial data: 



(3.41) 


t 


and one can read a from the asymptotics of the isomonodromic how near t = 1 [27, 30] : 


A = 1 + r/(f - 1)^-'"'^ + ... 


(3.42) 


valid when 0 < < 1. This particular asymptotic behavior at t = 1 will be clarihed 

in the next section. The constant r] is related to the other monodromy parameter 
Tr(MoMi) albeit by a complicated expression [2]. 

Asymptotics are more easily obtained for the r-function, dehned as: 



(3.43) 


is equal to the identity matrix at z = A. 
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which is related to the parameters of the dynamical system by 


^logr(t, {0i}) = + ^ + 


eA 

t - 1 


/s:i(A — t) 

tit-l) 


A(A — l)/i 


A(A-i)(A-t) r 2 
t{t - 1) f 


00 01 \ ^1/^2 
X^A-l^A-ty^^A(A-l) 


Poet OA 
t ^ t-V 
(3.44) 


where it is assumed that it'(t), A(t),/i(t) satisfy the equations of motion. Inspecting 
the parameters, we can arrive at the more direct correspondence [31]: 


A'(A(i).;«(«); t. So. 9.. 9,. 9o„) = 4 log T(i; 9„ 9, - 1.- 1) - ^4^:^ - 

(3.46) 

The r-function plays a central role in the theory of integrable systems, being interpreted 
in generic grounds as a generating functional, and its existence stems from a zero 
curvature condition. Despite the arguments, the r-function also depends on the trace 
of the composite monodromy operators M^Mt, MiMt and M^oMt. With the initial 
conditions set by (3.41), we have 


^(^- l)^logr(t,0,(T) 


t=to 


t=to 


toPtPi + (^0 ~ 1)000* + to(to ~ l)-^o 

(00 + 01 + Kl)0i = ~^{Po + 01 ~ 0* + 0oo), 


(3.46) 


where 0 = {0o , 01,0io,0oo} and a = {croi,ait} parametrize the invariant monodromy 
data - see Appendix A. The function ( = t(t — l)J^logr(t) obeys the second order 
differential equation: 


(i(t-i)C")' 




—2 det 


29(i 

tC'-C 


Vc' + 9S + 9? + e? 


t('-( C' + 9 ? + 9 ? + 9 ?- 9 ^\ 
29? (t - 1)C' - C 

el (i - 1)C' - C 29? / 

(3.47) 


sometimes called the “a-form” of Painleve VI equations [7, 31, 32]. Although the 
initial value problem is well-posed from the ODE perspective, with those conditions 
specifying an unique solution of the a-form of the Painleve VI equation (not to confuse 
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with the composite monodromy a), the second equation seems strange because it does 
not depend on the accessory parameter Kg. It seems to play the role of a consistency 
check on the hrst equation in (3.46). 

The conditions (3.46) are quite striking: from a formal point of view, the Hamil¬ 
tonian structure (3.40) allows for the complete solution of the monodromy problem 
using Hamilton-Jacobi techniques [33]. Due to the very peculiar integrable structure 
stemming from the Painleve property, however, the additional integration involved in 
hnding the generating function between A,/i and the canonical pair of monodromy 
variables [2, 34] is not necessary: it is already given by the r-function. 

3.2 The Painleve VI system near t = 0 and t = 1 

The r-function for the Painleve VI system is well studied. In a seminal work, its 
asymptotics were derived [30] , and, more recently, a proposal for the full series near the 
critical points was presented [6]. The asymptotics and the solution of the connection 
problem is enough information to completely dehne the r-function. We present the 
series near f = 0 and the asymptotics in the Appendix A. As seen above, given the 
parameters and the composite monodromies 

2cos7r(Toi = Tr(MoMt), 2 cosTraii = Tr(MiMt), (3.48) 

the r-function is uniquely determined. Then, formally, one could invert the relations 
(3.46) and hnd cr^ as functions of tg and Kg, as well as {Oi}. With this information, 
one can use the formula (3.35) to compute transmission coefficients. 

Approximate expressions for agt or au can be easily obtained to lowest order when 
to ~ 0 or to ~ 1; respectively. For the choice of coordinates (2.7), these limits corre¬ 
spond to the two near-extremal cases of Kerr-dS. Let us consider the to ~ 0 case hrst. 
The r-function^ given by [30] is 

r(t) oc ^ (3.49) 

where we assume that 0 < JJciot < 1, corresponding to the hrst terms of the expansion 

^Because of different definitions, our r function is related to Jimbo’s one by Tours(t) = — 

t)®i®*/^rjimbo(t)- 
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(A.l). Then we have that 


tit - 1)^ logr = - (^0 - Otf] + (1 - (rot)At^ + ..., (3.50a) 

to next-to-lowest order as t goes to zero. Applying the bonndary conditions (3.46), we 
get, to lowest order in to, 

The calculation for to ~ 1 is entirely analogous, using the expansion: 

r(t) oc (1 - + C>((1 - t)^^^“, (1 - t))], (3.52) 

yielding 

Cit = 6^1 + — 2 0 —^ (1 ~ ^o) + Ci((l ~ toY)- (3.53) 

As to —t 0, we are studying the “usual” extremal limit where the inner and outer 
horizons of the black hole coincide. Unfortunately, knowledge of aot only partially 
solves the problem in this regime, since the transmission coefficient (3.35) depends on 
the composite monodromies of the points involved in the scattering, in this case ^ = to 
and z = 1, corresponding to the outer and cosmological horizon respectively. In the 
Appendix A we calculate the relevant parameter in first non-trivial order in to. The 
second limit to —)■ 1 is related to the “large black hole” limit, where the outer horizon 
and the cosmological horizon coincide. Although the physical significance is not clear, 
the result is easily obtained and simple enough to worth the note. 

While the asymptotic expansion of au for small t can in principle be extracted from 
(A.l), there is an alternative way. The manifold of monodromy data is parametrized 
by seven numbers, 9 = {9q, 9^, 9i, 9ao} and a = {dot, ctoi, cxii}. However, because of the 
Fricke-Jimbo relation (see Appendix A), only six of those parameters are independent. 
Therefore, the manifold of “non-trivial monodromy data”, where the 9 are fixed, is 
two-dimensional. Moreover, this manifold is symplectic, the structure following directly 
from the Atiyah-Bott symplectic structure on the space of flat connections [26] - see 
[34] for the derivation of the formulae below and [2] for a review. A Darboux set of 
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coordinates a = a^t and 'ifj parametrizes the composite monodromies as follows: 


Pot =2 cos TTcr, 



(3.54) 


Pm = 7 ^. -7\/COiCloo - ;; (PoiPll - PoPl - PiPoo), 

2 sm Tiw 2 


sin TTcr ,_ 1 


where pt = Ti Mi, ptj = TrMiMj and 



On the other hand, the moduli space of flat connections has another set of Darboux 
coordinates stemming from the parameters of the Heun equation [12]: 


(3.56) 


= da A d'l/j = dK A dt, 


where O is a symplectic form. This then allow us to interpret the solution (3.46) in 
a different light: just as the derivative of the r-function with respect to t gives K as 
a function of a = aot, it is a solution of the Painleve VI Hamilton-Jacobi equation 
dehned by the Hamiltonian system (3.40). Therefore, the r-function is the generating 


functional of the canonical transformation between the two sets of Darboux coordinates 


(3.56). It follows then that 



which, in principle, gives a way of computing directly the monodromy parameters. 
However, the r-function dehned in (3.43) is dehned up to a constant which in principle 
could depend on the monodromy data. This means that the equation above (3.57) is 
dehned up to a function of 6 which could in principle be found from asymptotics, in 
a procedure similar to [33]. We will follow the more pedestrian approach of solving 
(3.46) in Appendix A. 
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4 Relation to Liouville conformal blocks 


The r-function (3.46) was solved combinatorially - see (A.l) - in the context of con¬ 
formal blocks in conformal held theories (CFTs) [6]. We will review the relationship 
between Fuchsian equations and correlation functions of primary operators in Liouville 
held theory. We will assume some familiarity with Conformal Field Theory, as in [35]. 
The subject is covered here as in [36], tacitly assuming the semiclassical limit. A more 
precise view can be seen in, for instance, [33] and references therein, although it should 
be said that in these treatments the r-function plays a secondary role. 

The Liouville held theory serves as a model of two-dimensional gravity, or of the 
scale mode of the metric in higher-dimensional Einstein-Hilbert Lagrangean. The Liou¬ 
ville mode in two dimensions is represented by a spin-zero held z), which nonethe¬ 
less has a classically anomalous transformation law: 


^'{z', z') = 0(^, z) + - log 
7 


dz' " 
dz ’ 


(4.1) 


which in turn requires a deformation of the stress-energy tensor 

T(z) = 


(4.2) 


~ and of its complex conjugate T{z)- for Q = The mode expansion of T{z) 

dehnes the Virasoro generators: 


T{z) = Y,L^z-^-\ L 


dz 

27ri 


T{z)z^+\ 


(4.3) 


which in turn satisfy the Virasoro algebra (we will momentarily not distinguish between 
Poisson brackets and commutators): 

[Lni Lm] (^ '^')Ln+m T f )^n+m) (^A) 

with c = 1 -|- 3Q^. Operators (functions) like Va(z,z) =: : - the :: represent 

normal ordering of the product of operators - are called primary operators because 

'^Quantum corrections make Q = 2/^ + We will assume the semiclassical approximation 7 —>■ 0 
throughout. 
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they transform nicely under conformal transformations: 






( 4 . 5 ) 


where we have the BPZ formula = a{Q — q)/2. We will omit the antiholo- 

morphic dependence from now on. This behavior is encoded in the Operator Product 
Expansion (OPE) between the primary operator and the stress-energy tensor: 

T{z)Va{w) ~ - - ^Va{w) +-, (4.6) 

where the tilde ~ means “up to regular terms”. From the dehnition of the Virasoro 
generators (4.3) one can see that the regular terms do not contribute to the action of 
the charges on the primary operators, a fact due to Cauchy’s theorem. 

Now let us consider the primary held V-^/ 2 {z). The state |x) = E_..),/2(0)|0) satishes 
the null state condition-, for an unitary representation of the Virasoro algebra where 
L\ = one has that the state 

(-^-1 + ^-^- 2 )| x ) ( 4 - 7 ) 

has zero norm. The requisition that (semi-classical) Liouville held theory is unitary 
forces this state to decouple from physical states. So, generically, 

+ 'ZL.2)V..,,^(z)X({Zi})) = 0 , ( 4 . 8 ) 

where X{{zi}) is a generic local operator. Using the OPE between T{z) and V_^/ 2 {z) 
one can see that 


L_YI_^I2{z) = d,V.,/2{z), L_2V.^I2 {z) =: T{z)V.^I2{z) : . (4.9) 


Now, if X{{zi}) is composed of four primary operators H* one can use the OPE 

between T{z) and each of the primaries to yield the Ward identity: 


d: + ^ 




hi 


2 ^ \ {z - ZiY z - Zi 


(r-x/2(.-) !!''«.(••)>= O' 


(4.10) 


2=1 
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where hi are the conformal weights of the Va-{zi) given, by the BPZ formula, and the 
accessory parameters Cj are given by: 


Ci = log(JJl4.(2:i)). (4.11) 

i=l 


One can then recognize the Heun equation (4.10) as the Ward identity for the 5-point 
correlation function involving V-^f 2 {z)- This correlation function, seen as a function of 
z, can be thought of as the “classical prohle” of the Liouville held in the guise 
in the presence of the primaries Vcn{zi). The accessory parameters of the Heun equation 
Ci are given in terms of the four-point function involving primary operators. Equation 
(4.11) is related to our solution (3.46) because our version of Heun equation (2.5) is 
slightly different, with a hrst order derivative term. As it turns out, the correlation 
function in (4.10) and the solution of (2.5) are related by multiplication of a simple 
function, a “s-homotopic transformation” like (2.6). Of course, the problem here as¬ 
sumes that the GET is unitary (and modular invariant) and so our application, with 
negative “conformal dimensions” hi = has to be taken as an analytical continua¬ 
tion of these dehnitions in terms of Liouville held [37]. It is an interesting open problem 
to see whether this theory makes sense on its own. 

The four-point function serves then as the generating function for the accessory 
parameters Cj. On the other hand, one can see that the structure of the four-point 
function stemming from Appendix A is compatible with the OPEs of two primaries 
being of the form 


nGZ 


J'(ai,a2,cr;n) 
{zi - Z2)'^ 


^cr+n (^ 2 )) 


(4.12) 


where (P^+n being an operator in the Verma module of the primary whose conformal 
weight is parametrized by a. From this expression we infer that a (modulo integer) 
has the interpretation of the intermediate channel of the “scattering process” of chiral 
vertex operators where (q;i,q; 2 ) is taken to ( 03 , 0 : 4 ). At least in unitary theories, the 
form of the functions should only depend on the representation theory of the Virasoro 
algebra. Indeed, the equation above can be seen as a Clebsch-Gordon decomposition 
of the tensor product of two Verma modules. The formulae exposed in the Appendix 
A were derived in the c = 1 case, and are believed to hold for generic parameters 
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of the Painleve VI r-function, and in this sense they are applicable to the black hole 
scattering problem, where the “conformal dimensions” are negative. It is not only quite 
impressive that this expansion also hold for the semiclassical calculation in Liouville 
field theory but also the most generic case provided by the Kerr-de Sitter black hole 
scattering. 

5 Discussion 

In this article, we have given expressions for the scattering coefficients of a conformally 
coupled scalar field in a four-dimensional Kerr-de Sitter black hole in terms of an 
implicit expression involving the Painleve VI r-function from (3.35) and (3.46). This is 
accomplished by using the same isomonodromy technique introduced in [3-5] to show 
the Painleve property of the isomonodromic deformation of a generic Fuchsian system 
[13]. By the same token, a careful choice of boundary conditions also imply that the 
eigenvalues of the angular equation are also given by the same Painleve VI r-function 
with coefficients appropriate to the problem (see (B.IO) in Appendix B). Moreover, the 
recent advances in the understanding of the Painleve VI r-function as c = 1 conformal 
block expansion, following from the proof of the AGT conjecture [14, 15], gives a 
combinatorial expression for this function [6, 7]. It is our belief that this completes the 
analytical resolution of the scattering coefficients. 

Pragmatically speaking, the consequences of this method for the study of superra¬ 
diance and quasi-normal modes are now within reach of the same numerical methods 
used to solve transcendental equations. The r-function expansion seems to systematize 
the recursion relations for the scattering coefficients based on “patching” methods [8]. 
Also, the Fuchsian form of the equations of motion for helds stems from the integrable 
(and singularity-free) solutions of the Einstein equations. So, there is a priori no rea¬ 
son why the methods outlined here would not work for generic spin perturbations, the 
Teukolsky master equation. The prospects for studies of electromagnetic and gravita¬ 
tional perturbations are particularly enticing due to their applications to astrophysics. 
We have already shown [16] that the same construction can be obtained in the zero 
cosmological constant limit, but now the relevant r-function is that of Painleve V. 

On a more abstract level, it seems startling that Virasoro representation theory 
methods are useful in studying Fuchsian-type of differential equations. The existence 
of a similar enhancing of the global conformal group SL(2, C) for extremal black holes 
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was noted and explored [38, 39] and arguments were given for the same happening 
for non-extremal black holes [40]. The picture arising from our results are however 
more intricate: one can use the isomonodromy flow to relate the scattering of helds 
at a particular non-zero value of the accessory parameter to - corresponding to the 
physical, rotating black hole - to a “confluent” Heun equation where to 0. In this 
limit, the position of the apparent singularity A also coalesces with a regular singular 
point (3.42) - which would correspond to an extremal black hole - but in general the 
accessory parameters diverge. The interpretation for this symmetry in terms of four- 
point functions is not clear, but it seems that, if there is a conformal description of 
the non-extremal black hole, the states and primaries involved will not be in the same 
SL(2, C) invariant state as the one used to describe the extremal black hole, as the 
accessory parameters diverge. We believe that the results presented here will not only 
be useful to further the studies in astrophysical applications but also will help clarify 
the more technical issues listed above. 
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A The Painleve VI r-function and asymptotics 

Here we present more information about Painleve VI r-function and the asymptotic 
expansion of au when to goes to zero. First, let us remind of the general expansion of 
the Painleve VI which has been given in [6, 7] in terms of the c = 1 conformal blocks. 
The expression is^ 

r(t) = aot + aot + n-,t), (A.l) 

nGZ 

^In the references, the monodromy parameters are defined with an extra factor of 2: (Jij}there 

I/2, (72^/2}here ■ 
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where the structure constants C are products of the Barnes functions (dehned by the 
functional relation G{z + 1) = V{z)G{z), r( 2 ;) being the Euler gamma): 


G{e,a) = 


rie e'=± + eV))G(l + 1(^1 + e6*oo + eV)) 


n6=±^(i+''^) 

and B have the structure of conformal blocks, given by the combinatorial series: 

A 


(A.2) 


(A.3) 


summing over pairs of Young tableaux A, p with 

j 

- 



Figure 2. A sample Young tableau A = {7,5,2,1} and the relevant quantities for the 
combinatorial r-function expansion. 




{{9t + a + 2{i — j))^ — 6 *q)((6*i + cr + 2(i — j))^ — 0^) 


n 


j)(A'. + Pi - i - j + 1 + cr)2 
{{Ot -a + 2{i-j)f - 9l){{0i - a + 2{i 


■3)f 


X 




16h2(i, j)(Ai + p'- - i - j + 1 - ay 


(A.4) 


where {i,j) denotes the box in the Young tableau A, Xi the number of boxes in row i, 
X'j the number of boxes in column j and hx{i,j) = Aj + A)- — i — j + 1 its hook length 
(see Fig. 2). 
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The parameter s in (A.l) is related to a special parametrization of the Fricke-Jimbo 
relation 

PotPitPoi + Pm + pit + Pm +pI+ Pt +pI+ PmPtPoo = 

{PoPt + PlPoc)P0t + {PlPt + P0Poo)Plt + {PoPl + PtPoo)P01 + 4, (A.5) 

where pt = TrMj = 2cos7r6'j and pij = Ti MiMj = 2 cos naij. This relation is valid for 
any group of four SL(2, C) matrices obeying the monodromy identity 


MoMtMiMoo = 1. (A.6) 

Following [29, 30], if we £x the 6i and aot the above relation can be parametrized in 


terms of s as 

{Pot ~ 4) Pit = Dt^+s + Dt-S ^ (A.7a) 

{pit - 4) Poi = Dy„+s + Du-s~^ + Dufl, (A.7b) 

with coefficients given by: 

Dt,o = Pot {poPi + PtPoo) - 2 {poPoo + PtPi ), (A.8a) 

Du,o = Pot {ptPi + PoPoo) - 2 {poPi + PtPoo ), (A.8b) 

Dt,± = 16 sin I iOt T crm + e6'o) sin | (^i =f am + , (A.8c) 

e=± 

Du,± = (A.8d) 

Solving the system (A.7) for s, we get 


s^(cos7r(6't =F am) — cos7r6*o)(cos7r(6*i =F aot) — costt^oo) 

= (cos 7r0t cos 7101 + cos 7700 cos 770^ ± i siu 77aot cos VrCToi) 

— (cos TT^o cos 7r6*i + cos 7r6*i cos TT^oo =F ^sinTTCTot cos7r(Tit)e^™°‘. (A. 9) 
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Finding in the Near-Extremal Case 

Now we have all ingredients to find au in the limit to 0. The first terms of the 
expansion (A.l) are 


r(t) ~ i-V4-(eo-0dV4 ^ C-S-h^-^ + Cot + sC+t^+'" + (A.IO) 

where a = aot and 

(Oj - 9? - a^nel - 0^, - a^) 

” 8cr2 

_ r^(l=P(T) -p-r r(l + |(e6*o + 6*t ± cT))r(l + |(e0oo + ± O')) 

^ r2(l±cr)^Al^r(l + i(e0o + 0iTc^))r(l + |(e0oo + 0iTc4))^- 

- iGtT(Tf){0lo - {diTcrf) 

16(T^(1 ± (t)^ 

Using this expansion in (3.46), we find to next-to-lowest order: 

t{t - 1)^ logr = -|[cr^ - {Oo - 9tY] + (1 - + ... , (A.12a) 

= —(1 — ^ ~ {9o — 9tY] + .... (A.12b) 


d 

dt 


i0-l)|logr 


Setting t = 
of a (3.51). 
known 


to above and applying (3.46), we get from (A.12a) the approximate value 
From (A. 12b), we can find the parameter s in terms of what is already 


s ~ 


(1-C4)^g- 

{9o + 6*1 + fi:i)6*t 


(A.13) 


then we just need to plug this result in (A.7) and this gives pu = 2cos7rcrii. 


B Angular Eigenvalues 


In Chambers-Moss coordinates, the Kerr-dS metric yields the following angular equa¬ 
tion in the conformally coupled case: 




du{P{u)duS{u)) + ( -2a u + Ce - ~ = 0’ (^-1) 
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with 

P{u) = {1 + au^){l — u^), Cl = +— 2maoj), (B.2) 

where a = a/L and = 1 + Our interest in this equation is primarily in the de¬ 

termination of the eigenvalues Ci. Approximate expressions for the Ci can be obtained 
using Fade (rational) approximants, as in [41], There, only the spin 2 case has been 
treated, but their method can be applied to any spin. In our case, the equation reduces 
to a Fuchsian equation with 5 regular singular points: 


6±i = =Fm, 6^i_ = ±i{ux^L — am) = ±0, 9^0 = 1- (B.3) 

a 

Upon the change of variables: 

X = , S'(x) = - a;oo)/(a;), (B.4) 

with to = (i + a)^/iia and Xoo = {i + a)/2i. The function /(x) satishes Heun’s equation 
in canonical form 

(Pf / 1 + 0 1-0 l+m \df / p+p- to(to - 1)^0 \ r m e-N 

dx^ \ X X —1 x — to/dx \x(x — 1) x(x — l)(x — to)/ 

and the accessory parameters given by 


Qo 


Aia 

X* 


p+ = 1 m, P- = 1 
[Ci + m + 1 a^ijn — 1) — 2acjx^) ■ 


From now on, we dehne 


(B.6) 


9o = 0, 9i = -0, 9^g = -m, 9oo = m. 


(B.7) 


The eigenvalue condition is that there exists a regular solution of the angular system 
at M = ±1 (corresponding to the north and south poles of the sphere), i.e. 


S{u) 


+ u ^-1 

{l-u)^/\ M^+1, 


(B.8) 
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which corresponds to the points x = and a: = oo in (B.5). This poses constraints 
in the connection matrices -E- 1,1 = Using the Wronskian normalization, (B.8) 

implies that one of the natural solutions at x = to (w = +1) is also a natural solution at 
X = 00 {u = —1). Thus, using (3.33), we conclude that the only allowed values for Ci 
will be those which the monodromy matrices and M^a are mutually diagonalizable. 
In fact, in this case, as the monodromies are integers, one expects logarithmic behaviour 
and the monodromies are only block diagonal. In any case, the composite monodromy 
<7*000 will satisfy: 

<7*000 = 2£, £ G Z, (B.9) 

and can in principle be obtained from (3.46) 

d 777 .© 

Qo{i) = ^logr(to; {0, -Q,m,m,atQ00 = 2(£ =f m), ao*^}) - (B.IO) 

as described in Appendix A. Thus the eigenvalue problem is also solved - somewhat 
formally - by isomonodromy. Asymptotic expansions for r near xq = 00 can be found 
in [7, 30]. 
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